Dynamical electroweak symmetry breaking due to strong Yukawa interactions 
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We present a new mechanism for electroweak symmetry breaking (EWSB) based on a strong 
Yukawa dynamics. We consider an SU(2)l X U(1)y gauge invariant model endowed with the 
usual Standard model fermion multiplets and with two massive scalar doublets. We show that, 
unlike in the Standard model, EWSB is possible even with vanishing vacuum expectation values 
of the scalars. Such EWSB is achieved dynamically by means of the (presumably strong) Yukawa 
couplings and manifests itself by the emergence of fermion and gauge boson masses and scalar 
mass-splittings, which are expressed in a closed form in terms of the fermion and scalar proper 
self-energies. The 'would-be' Nambu-Goldstone bosons are shown to be composites of both the 
fermions and the scalars. We demonstrate that the simplest version of the model is compatible with 
basic experimental constraints. 
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I. INTRODUCTION 



The majority of current models of electroweak sym- 
metry breaking (EWSB) can be divided into two large 
classes, according to the mechanism the symmetry break- 
ing is achieved (omitting exotic scenarios such as extra 
dimensions) : (i) The electroweak symmetry is broken at 
tree level by non-vanishing vacuum expectation valuc(s) 
(VEV(s)) of some elementary scalar ficld(s). This class 
includes the Standard model (SM) [l[ and its various ex- 
tensions like the two-Higgs-doublet models the min- 
imal supersymmetric Standard model or the little 
Higgs models [1, Q. (ii) The electroweak symmetry is 
broken non-perturbatively by strong dynamics of some 
new, yet unobserved, interactions. This class includes 
theories with new strong gauge interactions such as (ex- 
tended) technicolor [f| 0, HI, 0] , top condensate [ljj |Tl[ , 
or topcolor [H as well as effective models of the 
Nambu-Jona-Lasinio type [3, . 

In this paper we study a different scenario, which com- 
bines features of both classes mentioned above (elemen- 
tary scalars and strong dynamics): The electroweak sym- 
metry is broken non-perturbatively by strong Yukawa in 
teractions (a similar approach can be found in 16l 
and, in the context of non-relativistic systems, in 
GjlHil)- We believe that this option has certain advan- 
tages. For example, it allows to distinguish otherwise 
identical fermions at the level of Lagrangian by differ- 
ent Yukawa couplings and hence avoid dangerous (be- 
cause unobserved) Nambu-Goldstone (NG) bosons re- 
sulting from the spontaneous breaking of the large fla- 
vor global symmetries. The idea has already been stud- 
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ied in our previous paper [2l| using a prototype model 
with Abelian gauge symmetry. We have shown that the 
dynamically generated masses of different fermions may 
differ by orders of magnitude even with Yukawa coupling 
ratios of order one. We thus suppose that the present 
mechanism could shed some new light on the problem of 
the fermion mass hierarchy in SM. 

Encouraged by the previous model study, we investi- 
gate in this paper the possibility of EWSB by the same 
mechanism, considering for simplicity a single family of 
SM fermions. At this exploratory stage we do not want 
to make any kind of phenomenological predictions, our 
aim is merely to demonstrate that the extension of the 
scenario to the electroweak symmetry brings no new prin- 
cipial problems. In Sec. |TT1 we first introduce a model 
with global 5J7(2)l x U(\)y symmetry, which is bro- 
ken dynamically by the strong Yukawa couplings. In 
Sec. IIII1 we then gauge the symmetry, which induces dy- 
namical EWSB. In Sec. IIVI we present numerical results 
and perform a basic check of consistency with precision 
electroweak measurements. The outlook for future work 
as well as potential applications in other fields of physics 
arc discussed in the conclusions. 



II. THE MODEL 

We consider an SU(2)l x U(1)y invariant Lagrangian, 
equipped with one generation of the SM fermion multi- 
plets: The left-handed isospin doublets, denoted as = 
(vl,cl) T , <1l = {uL,d,L) T , and the right-handed singlets 
e R, URi d,R. In addition, we introduce the right-handed 
neutrino singlet isr in order to generate the neutrino 
mass. However, we will consider only the Dirac mass 
here, since the Majorana mass would require a more elab- 
orate formalism. We also introduce two complex scalar 
doublets, S = fi(°)) T and N = N^) T . Un- 

like in SM, they have ordinary masses Ms, M^ (i.e., 
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Mj N > 0), which implies they do not condense at tree 
level. Hence there is no need for the scalar self-couplings 
to stabilize the vacuum and they will for simplicity be 
neglected. 

Of key importance are the Yukawa couplings of the 
fermions to S, N, by assumption responsible for EWSB: 



£ Yukawa = Ve^L&R.S 

+ h.c. 



UuqLU R N 
(1) 



The Yukawa couplings can always be made real by a suit- 
able change of phases of the felds. For simplicity we do 
not consider interactions of the charge conjugated scalar 
doublets, S = \cr2S* and N = ic^iV*. Our 'minimal' 
Yukawa interaction Lagrangian in fact turns out to be 
sufficient for generating fcrmion and gauge boson masses. 

In the present paper we deliberately neglect renormal- 
ization and running of the Yukawa couplings. Within our 
model this is decently justified by the fact that thanks to 
the structure of the interaction Lagrangian ([I]), the one- 
loop pcrturbative corrections to the Yukawa couplings 
are absent. Throughout the calculation, the couplings 
are therefore implicitly assumed to be fixed and renor- 
malized at the scale of the scalar masses Ms,n- 

Of course, once the couplings become large and more- 
over the chiral symmetry is spontaneously broken, this 
naive argument breaks down. We do not attempt to de- 
termine the full nonperturbative flow of the couplings 
but merely note that it can be analyzed using the 
renormalization-group techniques, as was done for simi- 
lar models in Refs. [22|,[23j]. Our results, in particular the 
existence of a critical coupling for spontaneous symmetry 
breaking, are in qualitative agreement. 

With the above remarks in mind, we assume the full 
fcrmion propagators to have the form (/ = v, e, u, d) 



(2) 



where the S's are real functions of p 2 (the argument will 
be often omitted throughout the paper). Notice that the 
E's are scalar functions, which means that we also neglect 
the perturbative wave function renormalization. 

Likewise, for the full scalar propagators we adopt the 
Ansatz (X = S, N) 
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(®SN$SN 
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(3a) 
(3b) 



written using the Nambu doublet notation § x = 
(Y(°),Y(°)t)T $ SAr = (£(+),iV(-)t) T . The ITs are 

complex functions of p 2 . 

The spectrum can now be revealed by looking for 
the poles of the full propagators. In the fermionic 
case, provided there arc non- vanishing self-energies = 
T,fi + iSf 2) one arrives at the simple pole equations 
p 2 — |S/(p 2 )| 2 = 0. For the neutral scalar doublets <E>x 
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FIG. 1: Schematic diagrammatic representation of the 
Schwinger-Dyson equations (J5J. Full blobs represent the 
full propagators, the shaded blobs the one-particle-irreducible 
propagators proportional to the E's and II's, respectively. 
The letters L, R denote the fermion chiralities. 



the pole equation reads (p 2 ~ M x ) 2 - |IIx(p 2 )| 2 = 0. 
Its two solutions Mxi,2 are to be interpreted as that in- 
stead of one electrically neutral, yet complex scalar field 
Y(°) with the mass M x there are two real scalar fields 

J;] with distinct masses Mxi,2 in the spectrum. These 
mass eigenstates are apparently linear combinations of 
ReY^ ) and Iml' '. Similarly considering the charged 
doublet $sjv> its pole equation is [p 2 - M|)(p 2 - M^) — 
|ngjv(p 2 )| 2 =0. In this case the two solutions Msni,2 
indicate that the charged fields and N^~^ need to 
be mixed in order to get the mass eigenstates. In prin- 
ciple, charge conservation also allows mixing of S^ ' and 
iV' ' , but this only appears at the two-loop level. 

For later purposes it is convenient to define the 'anoma- 
lous' propagators, 



Sf( P ) 
Dx(p) 
Dsn{p) 



p 2 -|£/| 2 ' 
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IIsw 
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(4a) 
(4b) 
(4c) 



It should be emphasized that, once the symmetry is bro- 
ken dynamically by the above defined anomalous propa- 
gators, nothing any longer prevents the scalars from ac- 
quiring nonzero VEVs. In general, they will be induced 
by perturbative loop effects. However, since it is our aim 
here to demonstrate the possibility of dynamical EWSB, 
not driven by scalar VEVs, we will neglect these VEVs al- 
together in our Schwinger-Dyson (SD) equations. Their 
self-consistent incorporation will be studied in a future 
work. 

The self-energies S, n arc subject to the SD equations, 
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N = 



lis 



JV 



-2xY,Ncy 2 D J ^S D (k)S D (k + p) 
-2iY^N c y 2 u J ' -j^SuWSuik+p), (5d) 
-2iY,N c yuy D J ^Su(k)S D (k + p). (5c) 



(iVc denotes the number of colors). The sum in the last 
equation is over all isospin doublets, T> = (U, D) T . In our 
case, T> = £,q. The numerical solution of these equations 
is discussed in Sec. [IVJ 



III. GAUGE BOSONS 

Since the SU (2)l x U(1)y symmetry is broken down to 
electromagnetic E/(l) em , the question about the nature 
of the corresponding NG bosons (to become the longi- 
tudinal parts of the W^, Z bosons) arises. It turns out 
that, according to general principles, they can be 'seen' 
as massless poles in the proper vertices of the currents 
associated with the broken generators (24J. The residua 
of these poles are proportional to the symmetry-breaking 
proper self-energies £ and II. 

Without going too much into technical details (inter- 
ested reader can find them in [HI), ^ us only sketch the 
way leading to computation of the gauge boson masses. 
In this context it is useful to consider the vertex functions 
of the broken currents, 



Gijj a {x,y,z) 
G$ a (x,y,z) 



(0\T{j%(x)V(y)V(z)}\0) , (6a) 
{0\T{j£(x)*(v)&(z)}\0). (6b) 



Here j£ (a — 1,...,4) are the conserved currents as- 
sociated with the generators of SU(2)l x U(1)y and 
$ = $5, $jy, $5at. Taking the divergence of the vertex 
functions we arrive at the Ward-Takahashi (WT) iden- 
tities for the proper vertices of the gauge bosons with 
the fermions and scalars. For instance, for the proper 
vertex Zff we get (similarly for other gauge bosons and 
fermions / scalars) 



i(//V T /-W/"V : 



(7) 



FIG. 2: Contributions to the polarization tensor of a gauge 
boson (W or Z) from the fermion and scalar loops. The first 
vertex in each loop is bare, while the second one is the pole 
part of the proper vertex, extracted from the WT identity. 



with q — p'—p denoting the momentum of the incoming 
gauge boson and 



T f = 



2 cos Ow 



and Tf = joTj-fo- (T3 and Q denote the isospin and 
charge matrices, diagonal in the flavor space.) These 
WT identities remain valid even though the correspond- 
ing symmetries are spontaneously broken. Consequently, 
it is easy to see that upon inserting the full propagators 
©, © in the WT identities and taking the limit q -> 0, 
the proper vertices T M associated with the broken sym- 
metries (i.e., r^r ± , r£) develop a massless pole. It is 
to be interpreted as the propagator of an intermediate 
scalar excitation, bilinearly coupled to the gauge boson - 
the NG boson. In the case of Vt.j its pole part has the 

form (£/ = S/,i + i7 5 E/ j2 ) 



T 3/ (l-7 5 )-2Q/sin 2 



(8) 



r M - (V v) = — 

1 Z//pole^ >P> q 2 



■±f(p')T,+Tf± f (p) 



(9) 



Since the broken symmetry is gauged, the NG bosons 
decouple from the spectrum and manifest themselves 
only as the longitudinal polarizations of the W ± and Z 
bosons. By the Schwinger mechanism [25| . these poles 
in turn give rise to the gauge boson masses, which are 
determined by the corresponding residues [24| . Within 
our approximation, the polarization tensors are given by 
the graphs in Fig. [2] with the insertion of the pole parts 
of the proper vertices. 

As a net result, the gauge boson masses can be ex- 
pressed via the sum rules 
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Sc/| 2 + |Sd| 2 



(lid) 



Thanks to the fact that the dynamically generated self- 
energies are suppressed at high momentum, all integrals 
above are UV-finite. 



IV. RESULTS AND DISCUSSION 

The numerical analysis of the (Wick-rotated) SD equa- 
tions ([5]) shows that: (i) Non-trivial, UV-finite solutions 
do exist, (ii) The solutions are found only for relatively 
large values of the Yukawa coupling constants (of order 
of tens), (hi) Large ratios of fermionic masses can be 
accommodated while having the corresponding Yukawa 
coupling constants of the same order of magnitude. 

The point (iii) above is promising in the quest for real- 
istic fermion mass hierarchy. Because of the large param- 
eter space which needs to be scanned this has not been 
accomplished yet. However, some achievements, which 
suggest that it should be possible, have been made. First, 
we accommodated the hierarchy between the lepton and 
quark doublets. For y v = 63, y e w 84, y u = 65, yd = 90 
(and M| = 2, M N = 1) we found m v < m e = 0(1(T 4 ) 
and m u < md = 0(1O -2 ). (Note that all masses are 
expressed in the units of M/v-) Second, we managed 
to generate a large hierarchy within one doublet. Con- 
sidering only the leptons and 'neglecting' the quarks 
{Vu = Vd = 0), we found m e /m„ — 0(1O 2 ), calculated 
for y v w 50, y e = 80 (and again Mf = 2, M N = 1). 
Nevertheless, it should be emphasized that this lepton 
mass ratio will be significantly enhanced by the seesaw 
mechanism once the Majorana mass term is taken into 
account. 

The numerical results above were based on the exact 
solution of the Wick-rotated SD equations. In this nu- 
merical procedure one chooses some values of the input 
parameters (the Yukawa coupling constants and the bare 
scalar masses), solves the equations and checks whether 
the resulting spectrum (fermion and gauge boson masses) 
corresponds to the real values. It is, however, difficult to 
find in this manner the combination of the parameters, 
giving realistic spectrum; not only because of the huge- 
ness of the parameter space, but also since solving the 
SD equations even for a single set of parameters is rather 
time-consuming. 



Thus, in order to be able to invert the above procedure, 
i.e., to choose the spectrum first and then find the cor- 
responding parameters of the model, we relaxed the ex- 
actness of the solution and solved the SD equations only 
approximately. For that purpose we chose the Ansatz for 
the fermion and scalar self-energies in the form of step- 
functions with a common position of the step A (i.e., 
S/tp 2 ) nif #(A 2 — p 2 ) and U x (p 2 ) fi x 0(A 2 - p 2 )) 
and plugged them into the SD equations. This Ansatz is 
decently justified by the observation that the self-energies 
obtained by exact solution of the SD equations exhibit 
similar behavior: They are nearly constant in a wide 
range of momenta and then drop fast to zero around a 
scale common to all particles in the model. 

With an Ansatz on the shape of the self-energies, the 
SD equations can of course no longer be satisfied iden- 
tically. Because mass generation is a low-energy phe- 
nomenon, we demand instead that the left- and right- 
hand sides are equal at p 2 = 0. This results in a set of 
several algebraic equations for the unknown parameters 
of the model, which can be solved much more easily than 
the full integral equations ©■ 

According to the sum rules (fTUj) . the heavy quark fla- 
vors contribute most to the gauge boson masses. Like- 
wise, they also dominate the scalar self-energies (|5"c|l - 
([5e]) . which tie the individual fermion masses together. 
We therefore neglected for the moment the lepton dou- 
blet and treated the quark one as the t- and 6-quarks. 
We kept the gauge boson masses at the physical values. 
Due to numerical reasons, however, we could only fix the 
quark masses to the rather large values mt = 405 GcV, 
mi, = 10 GcV. Nevertheless, note that the ratio of these 
masses is still realistic. We then found the following so- 
lution: 



Hi 
Vb 
M N 
M s 
A 



375, 
235, 

433 TcV , 
283 TeV , 
71 TeV . 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 



This result shows again some important features. First, 
the Yukawa couplings are of the same order of magnitude, 
while keeping the realistic ratio of the quark masses in the 
doublet. Second, the scalar masses are much higher than 
the EW scale, which is gratifying, as explained above. 
Our numerical investigations suggest that if we managed 
to pull the quark masses down to the physical values, 
the scalar masses would even increase by a few orders of 
magnitude. 

The 'cut-off' A also deserves a remark. It corre- 
sponds to the energy scale at which the self-energies 
fall rapidly to zero. The fact that A <C Ms, Mm sug- 
gests that the scalar self-energies are practically vanish- 
ing for momenta corresponding to the scalar pole masses. 
Consequently, the scalar mass splitting is negligible and 
the iS-parametcr tends to be very close to zero. (For 
lis = IIjv = TIsjv = 0wc would have exactly 5 = 0.) On 
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the other hand, since A>m fl mj, the quark self-energies 
are almost constant (and non-zero) at the energy scale of 
fcrmion masses (under 1 TeV) . If we considered all three 
fermion generations together with a mixing, this would 
consequently lead to the CKM matrix being significantly 
close to the unitary form [26| . 



A. Compatibility with electroweak observables 

While the realistic fermion spectrum together with 
the Yukawa coupling constants can be presumably ac- 
commodated, it on the other hand brings the prob- 
lem how to keep the p-parametcr (defined as p = 
M^/(M§ cos 2 9 W )) close to 1. Note that if we set 
Vv = y e , Vu = Vd and Ms = M N , the SD equations 
(|5l) yield £„ = E e , S u = S d and IL S = = Usn 
and consequently p = 1 exactly. (Which is not sur- 
prising, since in this case the Lagrangian possesses the 
custodial symmetry.) Of course, in reality y v ^= y ei 
Vu 7^ Vd, so the custodial symmetry is broken by the 
fermions. However, for the values m„ < m e = 0(1O -4 ) 
and 7Ji u < TTid = 0(1O -2 ) mentioned above, we find 
p = 1.012. It is thus apparent that one can achieve re- 
alistic fermion mass ratios and still keep p reasonably 
close to one. A detailed analysis reveals that key role 
is here played by the scalars: If their bare masses and 
the self-energies are of the same order of magnitude (i.e., 
Ms ~ Mn and lis ~ II n ~ lis at) so that they approxi- 
mately conserve the custodial symmetry, they can render 
p close to 1, provided they are heavy enough so that they 
can overcome the effect of the custodial symmetry break- 
ing in the fermion sector. 

The fact that the scalars tend to be heavy is in fact 
reassuring. First, we do not have to deal with the usual 
hierarchy problem, that is, radiative instability of the 
masses of elementary scalars at the TeV scale. Of course, 
the problem is just postponed to a higher scale. In this 
sense we understand our model as an effective approach 
which is valid at energy scales around and below the 
heavy scalar masses. We do not pretend to have an ul- 
traviolet complete theory. 

Second, the new scalars must be heavy enough in 
order to avoid constraints from flavor-changing neutral 
currents (FCNC). Even though there are apparently no 
FCNC in the present model, for we only consider a sin- 
gle family at the moment, anticipating a future exten- 
sion to all three families of SM fermions, we can make 
at least a rough, order-of-magnitude estimate. With sev- 
eral fcrmion families there will be explicit flavor-changing 
Yukawa interactions, allowing, for instance, for the de- 
cay pr — > e^S^ -*. The virtual heavy scalar can sub- 
sequently decay as — ► e + e~. Our Yukawa interac- 
tions will therefore induce the flavor-chaning muon decay, 
pT — > e~e + e~, with the amplitude being roughly given 
by y 2 /Mg. (We assume that in the absence of fine tun- 
ing, all Yukawa couplings, including the flavor-changing 
ones, will be of the same order of magnitude.) The dom- 
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FIG. 3: The ^-parameter plotted for the special case Msi,2 = 
Mni,2 = Msni,2 = Mi, 2. Note that, according to the Parti- 
cle Data Group [l^j, S = -0.10 ± 0.10. 



inant muon decay channel, with branching ratio close to 
100%, is pT — * e~v e u ll ^ whose amplitude is analogously 
proportional to Gf- From here we infer the estimate 
RR{p~ -> e~e+e") ~ (y 2 /G F M^) 2 . Taking the current 
experimental limit [27|, BR(/i~ — > e~e + e~) < 10~ 12 , we 
find Ms/y > 10 2 5 TcV. Our numerical calculations re- 
ported above suggest that this constraint might impose 
some tension on our model, but it is certainly possible to 
satisfy. 

The very introduction of new scalars also affects the 
Peskin-Takeuchi ^-parameter (28[. In order to estimate 
the scalar contribution to it, we set for simplicity the 
scalar self-energies to be constant. The resulting S- 
paramcter can then be written as S = Ss + Sx + Ssn, 
where 
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(The p is just an arbitrary mass scale introduced for 
convenience; the total S'-paramctcr is independent of 
it.) Taking into account the previous discussion of 
the p-parameter and the scalar masses, we plotted the 
S'-parameter for the special case Msi,2 = Mni,2 = 
Msni,2 = -Mi, 2- The resulting S'-parameter (see Fig. [3]) 
meets the experimental bounds for any value of the 
mass ratio M1/M2 from 0.01 up to 100. When in this 
special case the ratio M1/M2 is far from one, the S'- 
parameter is well approximated by the simple formula 



S 



6tt ( 



hi 



Mi 

Mi 



On the other hand, for M 1 /M 2 



close to one the S-parameter behaves like 
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V. CONCLUSIONS AND OUTLOOK 

We have presented a new possible mechanism and 
the corresponding formalism for EWSB. We have intro- 
duced elementary scalars and argued that EWSB could 
be driven not by their VEVs, like in SM, but rather dy- 
namically, by their Yukawa couplings to fermions. If this 
was the case, then the 'would-be' NG bosons should be, 
as shown by general arguments, composites of both the 
fermions and the scalars. 

The presented model certainly has some advantages: 
(i) It distinguishes the fermion species already at the 
level of Lagrangian by their distinct Yukawa couplings 
and hence avoids the dangerous flavor symmetries, (ii) 
It ties the gauge boson masses up with the fermion and 
scalar masses by means of the sum rules (fTUj) . (iii) Since 
the resulting fermion masses are nonlinear functions of 
the Yukawa coupling constants, there is a possibility to 
accommodate large ratios of fermion masses, as seen in 
experiment, while keeping the Yukawa couplings of the 
same order of magnitude. This is supported by the nu- 
merical results shown above, as well as by the numerical 
analysis of our preceding model [2l| , dealing with a sim- 
pler toy version of the present model of strong Yukawa 
dynamics. 

So far we have considered only one fermionic genera- 



tion. Including all three generations will require to intro- 
duce the flavor mixing through the Yukawa coupling con- 
stant matrices. This opens a new way to study dynamical 
CP violation in both the quark and lepton sector, since 
even though the Yukawa coupling constants are real, the 
solutions to the SD equations are in general complex [2^ | , 
hence leading to complex 'mixing matrices'. 

Upon necessary modifications, the present mechanism 
has also the potential to describe non-relativistic bosonic 
superfluidity without a single-field condensate. The idea 
is supported by Refs. [l8[ an d [HI, in which the Cooper- 
type pairing in the many-boson system of 4 He is studied. 
The anomalous Green's function of the type ( ( I >( I > }, con- 
sidered here, was also studied in Ref. [2(|, along with the 
conventional Bose condensate. 
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